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1 Introduction 



Disentanglement is the process that transforms a state of two (or more) subsystems into 
an unentangled state (in general, a mixture of product states) such that the reduced 
density matrices of each of the subsystems are unaffaected. 

Let p*^"* be any entangled state of two qubits 1 and 2; and let pi, p2 be the reduced 
density matrices of 1 and 2 respectively. Then the operation of any disentangling machine 
(DM) is defined as 

^ent ^ pdisent 

together with 

p, = Tr,(p-*) = Tr.^p'^--*), i ^ j; i,j = 1, 2 

for all p^"*. This kind of ideal universal disentangling machine does not exist [Q. 

So the next question is whether there exists a disentangling machine which disentangles 
entangled state, and for which 

Tr,.(p^-«*) = ^^p^ + {^-^) /, z ^ j; ^, J = 1, 2 

where r/j (0 < r^j < 1 for i = 1, 2) is independent of p"^"* Recently it has been shown 
that m this kind of machine exists, by using local cloning operations, where the input 
states are all pure entangled states. Reference considered two cases, (1) r/i = 1 (or 
ri2 = 1), i.e., using only one local cloning machine, (2) r]i = 772 (= rj, say), i.e., which 
uses two local cloning machines with same fidelity. For the case (1), the maximum value 
of ri2 (or ?7i) is 1/3. In the case (2), the maximum attainable value of 77 is 1/V3. In the 
present paper, we want to find out the optimum values of rji and ri2, or, in other words, 
the optimal curve (if it exists) satisfied by t]i and ri2 {i.e., reduction factors corresponding 
to the optimal disentangling machine), by using most general (asymmetric) local opera- 
tions. Surprisingly, we got the same upper bounds on rj as has been found in in the 
corresponding cases (1) and (2). We have also obtained the optimal curve in the most 
general case, when asymmetric local operations are used. 

In section 2, for simplicity, we first consider disentanglement process by applying our 
disentangling machine on one of the subsystems. In section 3, we deal with the symmetric 
case where the same disentangling machine is used locally on the two subsystems. Next, 
in section 4, the most general disentangling machine is considered using asymmetric local 
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operations, where we discuss the disentanglement of mixed states. In section 5 we sum 
up our results and put some arguments regarding nonlocal operations. 



2 Totally asymmetric optimal universal disentangling 
machine 

In this section, we shall consider how we can disentangle a two qubit pure entangled state 
by local operation on any one qubit. Suppose we have two parties x and y sharing an 

entangled state of two qubits given by 

IV'> = «|00),^ + /3|11),^, (1) 

where a and (3 are non- negative numbers with + = 1. The first qubit belongs 
to X and the second belongs to y as usual. Now a universal transformation (a unitary 
operation) is applied to the qubit belonging to any one (say, y) of the two parties. This 
gives rise to a composite system PxyM consisting of the two qubits and a (disentangling) 
machine M . Tracing out on the machine states wc get a two qubit composite system 
which is disentangled {i.e., separable) under certain conditions. 

Consider the following unitary transformation Uj (associated with a machine state 
\M)j) applied on one subsystem j (where j — x or y), defined by 

Uj |0)^. |M)^. = mo, |0)^. |Mo)^. + my |1) . \M^). , (2) 

\M)^ = mo, |0)^. \Mo)^ + rfHj \M^), , (3) 

where ||Mo)j- , \Mi)j , Mq^ , Mi^ | are four normalized machine states, and (using uni- 
tarity) 

\mojf + \mijf ^ 1' 1 

Using orthogonality, we have (from (2) and (3)), 

mljffi^j.(Mo\Mo), + mljffHj.{Mi\Mi), = 0. (5) 

Now applying this operation Uj on an arbitrary one qubit state \4>)j = a\0)j + 
(where |ap + |6p = 1), we get the following composite state, 
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a [moj |0) . \Mo). + my |1) • |Mi)^.] + b [r^, |0) • |Mo) . + rr^j |Mi) 



amoj |0) . |Mo),. + amy |1) . |Mi) . + br^, |0), Mq) . + bfr^j |1) . Mi) . 



(6) 



We got (6) by applying the above unitary operation on any one qubit {x or y) pure 
state \(t>)j (j = x,y), where pj = \(p)jj{(t>\ = + (with \s\ = 1). And now we demand 
that the reduced density matrix, after tracing out the machine states in the equation (|[), 
is of the form |(1 + rjjS.a) (where < r/j < 1) for all s (isotropy) [Q. Then the machine 
has to satisfy the following equations : 



mo,mlj(M,\Mo). 



0, 
0, 
0. 



(7) 



Re m* mo, (Mo I Mo 



Re\mljm{j.(^Mi\Mi^^ 



0, 
0. 



r]j = mojniij .(Mi|Mo) . . 



(8) 
(9) 



\moj\ 

\rn^j\ 
\frHj\ 



2 



2 



1/2 



1/2 



(10) 



Now we apply the above-mentioned operation on one of the two qubits (say on y in the 
state \ip), in equation ([T|)0). Then the state = a |00)^^ + /5 \ll)^y is transformed (after 
tracing out the machine states, and applying all the above conditions i.e., (@), (^D, (0) ^ 
([T0|)) to the following density matrix : 



TA 



Q (l+'?y) n 
2 ^ 

2 





ia(3Ay 





-iaPAy 



'xyM\ — 

—iaPAy a[3rjy 

ia[3Ky 

/3'(i-^y) n 

2 ^ 

n /3'(i+%) 

^ 2 



^As the system x is unchanged, therefore, t]^ = I 
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where Aj = Im |mQjmoj-^.<^Mo|Moy .| (for j = x,y), and the entries of this matrix are 
arranged in accordance with the ordered basis {|00), |01), |10), |11)} of the two qubits. 

Now we apply the Peres-Horodecki theorem to test the inseperabihty of -Djjf , which 
states that a density matrix p of a bipartite system in (T ^ (T ^ is separable if and only if 
the partial transposition of p is positive semidefinite, i.e., each of the eigen values of the 
partial transposition of p is non-negative 0, 0, [0. It turns out that the state \ip) will 
be disentangled {i.e., D^^ is separable) if the following conditions are satisfied: 



universal disentangling machine), if ?7y < 1/3. 

Thus we see that all pure states of two qubits x and y will be disentangled by applying 
a disentangling machine locally on y, provided rjy < 1/3. 

Now that our requirement is also to have reduced density matrices d'"^] = Tiy[D'^^], 
^ad — Tr^l-Dj^^] of the disentangled state D^y as close as possible to those of the entangled 
one [i.e., pl^J = TTy[\tp){ip\] and pl^J = TTx[\ip){ip\] respectively), we note that the reduced 
density matrix of the subsystem x is unaltered whereas that of the subsystem y is changed. 
Let us summarise these results. 

1. It is possible to disentangle any arbitrary bipartite entangled state by applying 
local disentangling machine on one of its qubits provided the reduction factor (77) of the 
isotropic machine is less than or equal to 1/3. 

2. After disentanglement the reduced density matrices of the subsystems are given by 



where r]max = 1/3. 

3 Symmetric optimal universal disentangling machine 



«2/?2{(l+r/,)2(l-2r/,)-4A2}>0, (12) 
a'f3' {(1 - 3vy){l + Vyf + 8A^ (2A^ _ 1 + } > 0. _ 



All the three conditions in ([T^) will be satisfied for all a(3 G [0, 1/2] (as we are seeking for 





In the previous section we have shown how to disentangle any pure state of two qubits by 
applying local operation on one of the qubits. In this section we apply the local unitary 
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operation U = Ux = defined by equations @j and (|D, on both the parties x and y 
(in the state IV'), given in equation (P) separately. 

Each of the two parties now performs the same local unitary operation U on their 
own qubit, as described in the previous section. After this local operation, the reduced 
density matrix (tracing out the machine states) of the two parties x and y (applying all 
the constraints on the machine states, i.e., conditions (§), d^), (0) - (p^) is given by. 



Pxy 



+ a^T] - 2a/3A^ 
iaPAr] 
iaPAr] 



-ia(3Kri 









—iaPAr] 



iaPAr] 

-ia(3Ar] + p^T] - 2a(3A^ 



(13) 



It follows from the Peres-Horodecki theorem 0, 0, 0, that pxy is separable {i.e., the 
state \ip) is disentangled) if 



ai (2a2 + as) + as (aa + 2a3) - {a(3r]y{AA' + t]^) > 0, 

aia2 (a2 + 2a3) + ala-^ - (ai + as) {ai5r]f '(It^ + rf') 
- Aa2{a(3r]A)^ - 4{a(3fr]^A'^ > 0, 



(14) 



(15) 



so here i]^ = rjy ^ i] (say), = rriiy {i = 0, 1), = miy {i = 0, 1), \Mi)^ = \^'^i)y (* = 0, 1), 
A/A = mA = 0, 1), A, = A,;, = A (say) 



aia\a^ — 2a2 (ai + 03) (a/^r^A)^ — aia^i^aPYr]'^ — 2 (ai + 03) {aP)^ri^A'^ > 0, (16) 
where 
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(1 - ri)"^ 
ai = ^ — — ^ + a^r] - 2a/5A^ 



1 _ „2 
aa = + 2a/3A^ 

a3 = + _ 2a/5Al 

Now we shall consider the following two special cases, where in the first case, we put 
constraint on the machine, and in the second, we take the original state as a maximally 
entangled state. 

Case I: A = 0, i.e., Im{(Mo|Mo)} = [g. 

In this case, conditions ( p^ - ( jlB]) will be reduced to the following three conditions 
respectively: 

i^{3 + + 8(a/5r/)2} > 0, (17) 



(1 - r^y + 9,{al3ri)\l - 2r/^ - r]") > 0, (18) 

It is clear from the above conditions ( |T7|) - (|l^) that all bipartite pure entangled states 
{i.e., for all a/3 G [0, 1/2]), will be disentangled, if the reduction factor (r/) is less than or 
equal to 1/ Vs, and so the maximum value rjmax of rj is equal to 1/ "\/3-[] 
Case II:a(3= 1/2 g. 

Here, above conditions (p^ - (|1^) will be reduced to the following three conditions 



respectively: 



A'<^il-V'), (20) 

A' < (21) 

;i + + 4A2)(i + ^2 _ 4^2) _ 2^2 _ 3^4 _^ g^2^2 _ ^Q^4^ > Q ^22) 



As we have to maximize rj, we have to check (^) - (|2^) with all possible values of 
A. For that reason, we take, in the above three conditions ( [2(]| ) - (|2^), the values of A^, 



^Note the error in equation no. (17) of poj, and also in 
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starting from 0. Keeping all these in mind, it can be shown that all the conditions 
p^) will be satisfied for all aP G [0, 1/2], if the maximum value of r] is l/y/S. 



As mentioned in footnote 6, A = Im |mQmo <^Mo|Mo^| = — f]'^)/4: x 
Im{(Mo|iWo)} = ^/(l^rfjjl A, where A = Im{(Mo|Mo)}, so A G [-1,1]. In the 
general situation, for an arbitrary (universal) disentangling machine {i.e., for arbitrary 
value of A G [—1, 1]), we have to test whether the maximum value //max of the reduction 
factor can be made greater than 1 / -^/S. Now we note that the conditions given in (p!4D to 
(|T6|) are non-linear in A^, a (3 and rj, making it very difficult to get rjmax analytically from 
these conditions. And so we proceed numerically. As we are concerned with universal 
disentangling machines (different machines correspond to different values of A), therefore, 
we have to find out the maximum value rj^axi^'^) among all possible values ri{X'^) of the 
reduction factor for which all the states {i.e., for all the values of aj3 G [0, 1/2]) will be 
disentangled by a disentangling machine corresponding to the given value of A, so that 
all the conditions (|14D - (|16D are satisfied. Note that, our required rj^ax is the maximum 
value of all these r/max (A^)'s. From our numerical results, we have plotted ?7max(A^) against 
A^, in figure 1, which shows that the maximum value rjmax of all r/max (A^)'s occurs at A = 
(i.e., at A = 0), the maximum value being l/y/S. 



4 Asymmetric optimal universal disentangling ma- 
chine 

Here we apply the operations f4 and Uy (given in equations (|) and (|^)) separately on 
the qubits x and y of the state \ip) (given by equation (|l])) respectively. After taking trace 
over the machine states, and using the unitarity, orthogonality and isotropy conditions 
{i.e., equations @), (H), (0) - (0)) for both the parties x and y, the reduced density 
matrix (of the two parties x, y) becomes: 
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bi-c -ia(3A^r]y -ia(3AyT]^ al3rj^rjy 

iaPA^rjy 62 + c ia^Ayrj^ 

ia^AyT]^ 63 + c iaf^A^riy 

aPVxVy -i<y(3Ayr]^ -iaf3A^r]y - c 



(23) 



where 



b2 
b3 
b, 

c 



(1 - 7]^) (1 - riy) [t]^ + riy) 



4 „ 2 

{1 - T]^) {1 + r]y) ^ a' [T]^ - T]y) 



4 ^ 2 

(1-77^(1-^ ^ iVx + Vy) 
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(24) 



2aPA^Ay. 

Using Peres - Horodecki theorem [Q, 0, 0], we see that D^y will be separable {i.e., 
>) will be disentangled) if the following three conditions are satisfied: 



Fi + 2a(3 {r/.r/,A,,A, - a(3 (4A2 Aj + A^r^J _^ ^2^2^) | > 
F2 + (a/3)2 {AaPr],r]yA,Ay - [A^^ + Ajr/^) } > 0, 



(25) 
(26) 



{AAlAl (fi + r,l + SA^AJ + 2riir,f) + IGA^Aj (A^r^J + Ay^ ] 



a^(3^A^Ay7]^r]y 



[{2 -Til- ril - 16A2 A^) - 4 (Alril + A^^ } 



{4A^AJ [1+^1 + - ^tvf) + (AWy + K^l) (1 - Vlvl) + {AX - Kvl) {vl - n 
- \aPA,Ayi^,r^y (l - r^l) (l - r^j) > 0, (27) 



where 



F3 



6162 + 6163 + bibi + &2&3 + &2&4 + bsb4 - a^P'^rjlril, 
&1&2&3 + &1&2&4 + bib^bi + 62&3&4 - a'^P'^vlvl {h + b^) , 

6162&3&4 - bibia^l3'^rilrfy, 



(28) 



61, 62, ^3, ^4 being given by equations (p^). Here Aj = [1 ~ r]j) /A, where Aj 
Im(,(Mo|Mo).j for J =x,2/. 
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Let us first consider the case when = 0, and A„ = 0. In this case, above three 



conditions (|25|) - ( pT]) will be reduced to the following three conditions respectively. 

Fi > 0, (29) 

F2 > 0, (30) 

Fs > 0, (31) 

where Fj's are given in (^). Conditions ( PUD - (^) will be satisfied for all aP G [0, 1/2] 
if the reduction factors rj^ and rjy satisfy the following relation: 

1 , , 

VxVy < g- (32) 

Thus for given any i]^ {%) in (0,1], the maximum value r]y(max) {Vx{max)) of % (77^.) is 
1/37]^ (1/3%). 

Next we consider the general situation. Since we are looking for universal disentangling 
machine(s), thus, for given any value of the pair {Xx,Xy) in [0,1] x [0,1], and for given 
any value of rj^ G (0,1], each of the three conditions ( ^5]) - (pTl) has to be satisfied 
for all a/5 G [0,1/2], and for all values of rjy G (0, //^(max)], where < riy{max) < 1- 
Now we see that (i) satisfaction of the condition (^) (universally) implies satisfaction 
of condition (|29|) (universally) (as the term other than Fi on the left hand side of (pSf ) 
becomes non-positive for ajS = 1/2), (ii) satisfaction of the condition ( [26| ) (universally) 
implies satisfaction of condition (^) (universally) (as the term other than F2 on the left 
hand side of (^) becomes non-positive for a(3 = 1/2). But it can be shown that one 
cannot find even a single value of a/5 G (0, 1/2] for which the term on the left hand side 
of (P7|), other than F3, always remains non-positive for every choice of A^;, Xy G [0, 1], and 



for every choice of fjxjVy ^ [0? l]-0 So we have to take into account the three conditions 
( |25| ) - ( |27| ) all together. We have obtained numerically that for given any value of the 
pair (A^:, A^) in [0, 1] x [0, 1], and for given any value of rj^ G (0, 1], the maximum value 
Vy(max) (say) of rjy G (0, 1], for which all the three conditions (^Sf) - (^71) are satisfied for 
all a/5 G [0,1/2], satisfies the relation rjyi^max) < l/^^x (see figure 2 for a comparison). 
So the optimal disentangling machine, in the asymmetric case, will be obtained from the 
^We have verified this numerically. 
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case where A^^ = and = 0, and so the optimal curve to be satisfied by the reduction 
factors T]x and rjy, is given by the following rectangular hyperbola. 

VxVy = ^- (33) 

It is clear from condition (|3^) that (i) the maximum value of rjy, in the totally asymmetric 
case [i.e., when rjx = 1) is 1/3 (as shown in section 2), and (ii) the maximum value of rj, 
in the symmetric case {i.e., when rj^ = rjy = i]) is (as shown in section 3). 

Now we come to the question of disentanglement (by using local operations only) of an 
arbitrary mixed state p of the two qubits x and y. In its spectral representation, p takes 
the following form: 

p = j2p,pm], (34) 

i 

where Pi > 0, J2i f^i = ^, and P [I'j/'j)] is the projection operator on the (normalized) pure 
state \ipi) of the two qubits x and y. One can express \ipi), in its Schmidt form, as 

m=0^im),y+k\hh),y, (35) 

where a,, bi are non-negative numbers with af + hf = 1, and |li)^- are two orthonormal 
states in the Hilbert space of the qubit j (for j = x,y). As discussed earlier, each of the 
states will be disentangled by using local operations, with some common values 
of the reduction factors rj^ and rjy. Let p'^f^^) and p^ be the single particle reduced 
density matrices of \ipi) and p respectively, corresponding to the qubit j (j = x,y), before 
disentanglement; and let p'^f^-^ and p"*^ be the single particle reduced density matrices of 
\ipi) and p respectively, corresponding to the qubit j (j = x,y), after disentanglement. If 
r]j is the reduction factor associated with the disentanglement of p, corresponding to the 
qubit j (j = x,y), we then have 

Pf - + = E /^.P,%,) + (36) 

or 

E P^Pm) - E /^^ [v^Pm) + = ? + ^^^^ 

Thus we get from equation (|37|) that rjj = rjj for j = x, y. 
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5 Discussion 



Since an ideal universal disentangling machine does not exist, we have explored how well 
a universal disentangling machine can be, using most general type of local operations. 
We have shown that in the case of optimal universal disentangling machines, reduction 
factors lie on a rectangular hyperbola. Here one should mention that in the case of optimal 
universal clonning machines, the relation between the reduction factors follows directly 



from no-signalling phenomenon [[T^. So one may look for the physical phenomenon (or 
phenomena), from which above-mentioned rectangular hyperbola would follow directly. 
This is still an unresolved problem. 

We now address the issue of obtaining a better disentangling machine, if possible, using 
non-local operations. 

For all (non-unitary) local operations, entanglement decreases, whereas, for non-local 
operations, it may increase, decrease, or remain same. In the disentanglement process 
described here, a separable state remains separable without any further constraints (other 
than isotropy condition) on the machine, but even in order to keep a separable state 
separable under non-local operations, some constraints other than isotropy condition have 
to be imposed on the machine, which seems to decrease the reduction factor. But if one 
redefine the notion of disentanglement as a process which also keeps every pure bipartite 
product state, a product of two (single particle) density matrices of the two particles, 
then the allowed class of disentangling machines will comprise of only local operations. 

In this regard, for intuitive understanding, we point out that local cloning machine, with 
the blank copy, functions as universal disentangling machine which can be made optimal. 
Now if we use optimal local cloning machine, which produces three copies instead of two, 
then the cloning machine along with the two blanck copies acts as a universal disentangling 



machine with reduction factor being 5/9 |T^. A non-local cloning machine ||10[ along with 
six blanck copies, which produces seven copies of the bipartite states, acts as universal 
disentangling machine with reduction factor being 11/35 |]T^, which is much less than the 
former case. Recently Mor and Terno [|l^ obtained a set of bipartite entangled states, 
which can be perfectly disentangled, provided the reduced density matrices of one of the 
parties commute. But this is achieved using local operation, namely local broadcasting. 
In conclusion, we have obtained the optimal disentangling machine exploiting the most 
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general local operations, and discussed that there can not be any non-local operation 
which may give a better one. 

Acknowledgement: The authors thankfully acknowledge the sugessions and comments 
of the anonymous refree for revision of the earlier version of the manuscript. 



13 



References 



1. T. Mor, Phys. Rev. Lett. 83 No. 7, (1999) 1451. 

2. D. Terno, Phys. Rev. A 59 (1999) 3320. 

3. D. Brufi, D. P. DiVincenzo, A. Ekert, C. A. Fuchs, C. Macchiavello and J. A. Smolin, 
Phys. Rev. A 57 (1998) 2368. 

4. S. Bandyopadhyay, G. Kar and A. Roy, Phys. Lett. A 258 (1999) 205. Lett. A). 

5. A. Peres, Phys. Rev. Lett. 77 (1996) 1413. 

6. M. Horodecki, P. Horodecki and R. Horodecki, Phys. Lett. A 223 (1996) 1. 

7. M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 80 (1998) 2261. 

8. Here we have taken the identical values mo^: = mQy as mo, which can be taken as the 



real quantity J + v)/'^ (see equation (0)); similarly we have taken the identical 



values moa; = ruoy as tuq, which can be taken as the real quantity y{l— rj)/"^ (see 
equation ([10|)). Also we denote here the identical states \Mq)^ = |Mo) as \Mq), and 



the identical states 



Mr 



= Mo)^ as |Mo). Thus A = Im {mj^m^, (Mo|Mo)} = 
implies that Im |<^Mo|Mo^| = 0, as our requirement is to obtain universal disen- 
tangling machines, and it is known that such machines may exist provided t] is less 
than 1 [|. 

9. In this case, the reduced density matrices p'j'^ and pj'^ (of the maximally entangled 
state :^(|00)a;y + \ll)xy), before and after disentanglement respectively) are both 
equal to (1/2)/, / being the 2x2 identity matrix. So there is no direct effect of the 
reduction factor t] on the disentanglement of -^{\00)xy + |ll)xy)- But we try to see 
here how the universal disentangling machines work on this maximally entangled 
state. 

10. V. Buzek and M. Hillery, Phys. Rev. Lett. 81 (1998) 5003. 

11. S. Ghosh, G. Kar and A. Roy, |quant-ph/9907001| (to be published in Phys. Lett. 
A). 

14 



12. S. Bandyopadhyay and G. Kar, Phys. Rev. A 60 (1999) 3296. 

13. T. Mor and D. Terno, |quant-ph/ 990703^ . 



FIGURE CAPTION : Figure 1 describes the decreasing nature of y]max{^^) with the 
increament of A^. This graph shows that the maximum value rimax of Vmaxi.^^) occurs at 
= 0, and (so) rf^ax is equal to 1/v^. 

FIGURE CAPTION : Figure 2 shows the optimal curve rj^riy = 1/3, which corresponds 
to the case A^, = 0, Ay = (represented in the figure by the continuous line). Also, 
for a comparison, the optimal curves corresponding to the cases (A^,, Xy) = (0.2, —0.2) 
(represented in the figure by the broken line of the type '- - - -'), (A^., A^) = (0.5, —0.5) 
(represented in the figure by the broken line of the type -'), {\x,\y) = (0.9,0.1) 
(represented in the figure by the broken line of the type ' ') are given here. 
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